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Abstract 

Polarizabilities, dispersion coefficients, and long-range atom-surface interac- 
tion potentials are calculated for the n = 2 triplet and singlet states of helium 

using highly accurate, variationally determined, wave functions. 
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The advent of doubled basis sets has made it possible to calculate precisely many proper- 
ties of two-electron atomic systems . We apply variational methods developed previously 
and demonstrated for the helium atom ^ to calculate nonrelativistic values of the electric 
dipole, quadrupole, and octupole polarizabilities and corresponding dispersion coefficients 
for the metastable n = 2 singlet and triplet states, respectively, He(2^S') and He(2'^S'). Ad- 
ditionally, potentials for the atom-wall interaction of a He (2 ^5) or a He(2^S') atom and a 
single perfectly conducting wall or a dielectric wall are calculated with the inclusion of re- 
tardation effects due to the finite speed of light. Our results for atom- wall interactions are 
germane to experiments involving atom-evanescent wave mirrors 0. 

In this paper the notation of Ref. [0 is followed very closely; references to equations of 
Ref. [0 will be preceded by the symbol I. Atomic units are used throughout. 

The dispersion interaction of two like atoms can be written 

U{R) = -CeR-^ - CsR'^ - CioR-^°, (1) 

where the coefficients Cq, Cg, and Ciq are the van der Waals coefficients, R is the interatomic 
distance, and 

C6 = (3/7r)G(l,l), (2) 
C8 = (15/7r)G(l,2), (3) 
Co = (28/7r)G(l, 3) + (35/7r)G(2, 2) , (4) 

with 

POO 

G{l,m) = / ai{iuo)am{ioj)duo , (5) 

where aiiioj) is the 2'-pole dynamic polarizability function evaluated at imaginary frequency 
defined by Eqs. (6)-(9) of Ref. 0, and similarly for am{i^jj)- 

When the effects of retardation due to the finite speed of light are considered the potential 
U{R), Eq. (1), can be replaced by gi 
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V{R) = -CMR)R-^ - Csfs{R)R-' - Cio/io(i?)i?-'", (6) 

The coefficient /lo(-R) will not be considered in this paper as the Ciq term is usually neg- 
ligible. Expressions for the retardation coefficients, feiR) and fs{R), as integrals involving 
the dynamic electric dipole polarizabilities, are given in Eqs. I-(5) and I-(7). 

The form intrinsically includes certain relativistic effects, so that when fe{R) and 
fsiR) are expanded in powers of the fine structure constant afs = 1/137.035 989 5 for small 
distances 

V{R) ~ -i?"^[C6 - alR^Wi] - R-^[Cs - alR^We], (7) 



where 



Wi = - I duj uj^aliiu) (8) 
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and 
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Wq = — duj uP'ai{iLo)a2{iLo)- (9) 
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The relativistic origin of the coefficient W4 has been discussed by Power and Zienau see 
also |]11| . The coefficient of the factor q;|/ R^ in corresponds to the theory of Power 



and Thirunamachandran H and is equal to the coefficient Wll;a,2 in the theory of Meath and 



Hirschfelder based on the Breit-Pauli Hamiltonian |TT|. As the distance increases retardation 
arising from the finite speed of light becomes important and the potential approaches its 
asymptotic form, see Eqs. I-(13) and I-(14), 

V{R) ~ -KjR-^ - KgR'\ (10) 

with 

K, = = 250.81 K, = ^^iW^O) ^ ^^^^ ^ «,(o)a,(0). (11) 

An expression for the potential Vmd{Ri^) for the interaction of an atom and a 

dielectric wall was presented in Eq. I-(15), where R is the atom-wall distance and e is the 



dielectric constant of the wall. The expression is a double integral that can be evaluated 
with knowledge of the function ai{iuj). For small distances VAtoiR, e) has the limiting form 

^At.(i?,e)~-§^, (12) 

where 

C3 = — / dujai{iuj). (13) 
4:71 Jo 

As the separation increases retardation becomes important and the potential approaches its 
asymptotic form, 

l^AtD(i?,6)~-|^^0(6), (14) 

where 0(e) is given in Eq. I-(21) and 

= 3ai(0)/(87rafs) = 16.357ai(0). (15) 

For a perfectly conducting wall VAtoiR,^) reduces to VAtuiR), where 

VAtMiR) = VAtDiR, 00) = ~C^h{R)R-^ (16) 

and the retardation coefficient fz{R) is an integral involving ai{iuj) and is given in Eq. I-(26). 
For small distances VAtM(-R) ~ —C-i/R? and for asymptotically large distances VAtAf(-R) ~ 
—K^^/R^. Table 11 of Ref. |]14| summarizes the various limits of VAtoiR-i e)- 



It has been shown that double basis sets work well for calculations involving S states of 
helium 0. The basis set used here was constructed as in Ref. with basis set functions 
expressed using Hylleraas coordinates 

{x^Ao^.P) = r\TiT\,e-^^^-P^-}. (17) 

The explicit form for the wave function is 

*(ri, Y2) = Xijkiai, A) + 4ff, Xijk{a2, ^2)] ± exchange , (18) 

ijk 



and i + j + k < Q. The convergence of the eigenvalues is studied as Q is progressively 
enlarged. Finally, a complete optimization is performed with respect to the two sets of 
nonlinear parameters ai, jSi, and 0:2, (32 by first calculating the derivatives analytically in 



H 



oa / \ oa / 



(19) 



where a represents any nonlinear parameter, E is the trial energy, H is the Hamiltonian, and 
= 1 is assumed, and then locating the zeros of the derivatives by Newton's method. 
These techniques yield much improved convergence relative to single basis set calculations. 
The method of the evaluation of the two-electron integrals in Hylleraas coordinates can be 



found in Ref. 15 



The expressions for the dynamic dipole polarizabilities, Eqs. (6)-(9) of Ref. 0, were 
evaluated using the wave functions determined by the variational method. Values of the 
static polarizabilities are given in Table | for He(2 ^S") and B.e{2^S). The polarizabilities 
given in Table | are extrapolated results, with the convergence studied as in Refs. and 
and the estimated extrapolation error in the last digit is given in parentheses with 
the listed values. The largest basis set sizes used consisted of 616 functions for the S 
states, 910 functions for the P states, 931 functions for the D states, and 1092 functions 
for the F states. The converged results are compared with some previous calculations and 
experiments in Table || and|Tl[ Ekstrom et al. 0] determined the He(2 ^S) polarizability by 
combining their measured Na polarizability with the Molof et al. |1^ measurement of the Na 
polarizability relative to the polarizability of He(2 ^S). For the triplet state the experimental 
values of Ref. and of Refs. fl^Jl^l and the bounds of Glover and Weinhold are compared 



with our calculated polarizability in Fig. ||. 

The dynamic polarizability functions were constructed using the largest basis sets of 
each symmetry and used to evaluate the atom-atom dispersion constants and retardation 
coefficients. Our results for the dispersion constants are given in Table |IV|, with the estimated 
convergence errors given in parentheses, and the results are compared to other calculations 
in Tables and 0. The retardation coefficients are given in Table ^ and Fig. |. 



Chen and Chung [|T^] calculated the coefficients W4 and Wq for He(2 ^5*) and their results 
are compared with ours in Table |V111| ; their published value of Wq was multiplied by the 
factor I to correspond to the theory of Ref. |^ and the agreement is very good. 

For the atom-wall interactions the values of the coefficients C3 can be obtained from the 
alternate expression 

1 



/ N 
u=l 



, (20) 



which follows from integration of Eq. (]T3|) , where is the number of electrons and |0) is 
accordingly the 2 ^5* or the 2 '^S wave function. Since high-precision matrix elements are 
available i,|20| Eq. (^) was used to obtain the coefficients ^3(2 ^S) = 2.671 212 717025 and 
^3(2 = 1.900 924 084 097. 

The dynamic dipole polarizability was used to evaluate the potential for various dielectric 
walls. Results for He (2 ^S) are given in Table |I^ and Fig. ^ and those for He (2 ^S) are given 
in Table ^ and Fig. ^. The dielectric materials represented in the tables correspond to 
fused sihca (e = 2.123), BK-7 glass (e = 2.295), and a GaAs-type material (e = 3.493). 
The tabulated potentials may be helpful in planning and analyzing experiments with atom- 
evanescent wave mirrors, see for example Ref. 

We thank Professor G. W. F. Drake and Dr. P. L. Bouyer for helpful communications. 
The Institute for Theoretical Atomic and Molecular Physics is supported by a grant from the 
National Science Foundation to the Smithsonian Institution and Harvard University. ZCY 
was also supported by the Natural Sciences and Engineering Research Council of Canada. 
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TABLES 

TABLE L Values of the static polarizabilities ai(0), 02(0), and 03(0) for the 2^5 and 2^5* 
states of He. Numbers in parentheses represent the estimated error in the last digit of the listed, 
extrapolated value. 

State ai(0) 02(0) 03(0) 

2^5' 800.316 33(7) 7106.053 7(5) 293 703.50(6) 

235" 315.631468(12) 2 707.877 3(3) 88 377.325 3(7) 



TABLE IL Comparison of static multipole polarizabilities ai(0), 02(0), and 03(0) for He(2 ^S). 
For the experimental value numbers in parenthesis give the quoted error. 



Author (year) 



Ref. 



ai(0) 



02(0) 



03(0) 



Crosby and Zorn (77) Expt. 
Chung and Hurst (66) 
Drake (72) 
Chung (77) 

Glover and Weinhold (77) 
Lamm and Szabo (80), EGA 
Rerat et al. (93) 
Chen (95) 
This work 



ii 

m 



729(88) 
801.9 
800.2 
801.10 

803.31±6.6P 
790.8 
803.25 
800.34 

800.316 33(7) 



6870.9 



7106.053 7(5) 293 703.50(6) 



'^Bounded theoretical value. 
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TABLE III. Comparison of static multipole polarizabilities ai(0), 0:2(0), and 03(0) for He(2^S'). 



Author (year) 



Ref. 



ai(0) 



Q2(0) 



«3(0) 



Crosby and Zorn (77) Expt. 

Ekstrom et al. (95) Expt. 

Bishop and Pipin (93) 

Rerat et al. (93) 

Glover and Weinhold (77) 

Drake (72) 

Chung (77) 

Chen and Chung (96), B Sphne 
Chen and Chung (96), Slater 
Chung and Hurst (66) 
Chen (95) 
This work 



m 

Ml 

m 



122 
1 



301(20) 
322(6.8) 
315.631 
315.92 

316.24±0.78^ 

315.608 

315.63 

315.63 

315.611 

315.63 

315.633 

315.631468(12) 



2 707.85 
2 662.02 



2 707.89 
2 707.81 



1377.2 



1377.4 
1356.2 



2 707.877 3(3) 88 377.325 3(7) 



^Bounded theoretical value. 



TABLE IV. Values of Cg, Cg, and Cio for the interaction of two He atoms. 



System 


Ca 


Cs 


Cio 


2 


11241.052(5) 


817 250.5(4) 


108167630(54) 




3 276.680 0(3) 


210 566.55(6) 


21786 760(5) 
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TABLE V. Comparison of Cg, 


and Cio for the He(2 ^S')-He(2 ^5") system. 


Author (year) 


Ref 




Ce 




Cio 


Glover and Weinhold (77) 


M 




11330±630^ 






Rerat et al. (93) 


m 




11360 


812 500 




Victor et al. (68) 






11 300 






Lamm and Szabo (80), EGA 


m 




10 980 






Ghen (95) 


m 




11244 


817 360 


108184000 


This work 






11241.052(5) 


817 250.5(4) 


108167630(54) 


^Bounded theoretical value. 












TABLE VL Gomparison of Cs 


, and Cio for the He(2 35)-He(2 ^5) system. 


Author (year) 


Ref. 


Ce 




Cio 


Glover and Weinhold (77) 


[29| 


3 289±90'' 






Victor et al. (68) 


1 


m 


3 290 






Lamm and Szabo (80), EGA 


1 


m 


3 300 






Rerat et al. (93) 


1 


m 


3 279 


208 600 




Bishop and Pipin (93) 


1281 


3 276.677 


210 563.99 


21 786 484 


Ghen (95) 


[31| 


3 276.1 


210 520 


21 783 800 


Ghen and Ghung (96), B spline 


1 


19| 


3 276.10 


210 518 


21 783 800 


Ghen and Ghung (96), Slater 


1 


m 


3 275.90 


210 507 


21 780 200 


This work 






3 276.680 0(3) 


210 566.55(6) 


21786 760(5) 


^Bounded theoretical value. 
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TABLE VII. The dimensionless retardation coefficients fe^R) and fsiR) for tlie atom-atom 



interaction. The dispersion coefficients 


Cg and Cs are 


also g 


;iven. In the last line. 


labeled "Asymp- 


totic" 


the values calculated using the asymptotic forms /g 


- Kr/iRCe) and /g 


~ Kg/{RCs) are 


given 


in, respectively cols. 2,4 and cols 


. 3,5, with 


and Kg given in Eq. (11). 






He(2 i5)-He(2 ^S) 




He(2 35)-He(2 35) 






Cs 










112 41.052(5) 


817 250.5(4) 




3 276.680 0(3) 


210 566.55(6) 


R 


h{R) 


fsiR) 




h{R) 


fsiR) 


10 


0.999998 


0.999996 




0.999995 


0.999992 


15 


0.999996 


0.999992 




0.999988 


0.999982 


20 


0.999993 


0.999986 




0.999980 


0.999969 


25 


0.999989 


0.999978 




0.999968 


0.999952 


30 


0.999984 


0.999968 




0.999955 


0.999931 


50 


0.999958 


0.999913 




0.999879 


0.999812 


70 


0.999919 


0.999833 




0.999770 


0.999638 


100 


0.999840 


0.999666 




0.999548 


0.999281 


150 


0.999655 


0.999271 




0.999034 


0.998446 


200 


0.999408 


0.998742 




0.998358 


0.997337 


250 


0.999106 


0.998088 




0.997533 


0.995980 


300 


0.998750 


0.997318 




0.996573 


0.994397 


500 


0.996857 


0.993223 




0.991568 


0.986162 


700 


0.994325 


0.987790 




0.985052 


0.975560 


1000 


0.989563 


0.977772 




0.973189 


0.956657 


1500 


0.979675 


0.957749 




0.949690 


0.920633 


2000 


0.968032 


0.935332 




0.923467 


0.882350 


2500 


0.955170 


0.911784 




0.895919 


0.843974 


3000 


0.941459 


0.887865 




0.867915 


0.806627 
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5000 


0.882570 


0.795572 


0.759993 


0.675026 


7000 


0.822962 


0.714364 


0.666219 


0.572701 


10000 


0.739435 


0.614288 


0.554000 


0.460880 


15000 


0.622323 


0.492124 


0.424424 


0.342480 


20000 


0.530963 


0.407029 


0.340027 


0.270048 


25000 


0.459732 


0.345261 


0.282000 


0.221929 


30000 


0.403520 


0.298792 


0.240136 


0.187921 


50000 


0.266133 


0.191769 


0.149197 


0.115667 


70000 


0.196376 


0.140118 


0.107696 


0.083251 


100000 


0.140107 


0.099381 


0.075824 


0.058519 



Asymptotic 

100000 0.142912 0.100738 0.076257 0.058760 



TABLE VIII. The coefficients W4 and Wq appearing in the expansion of the atom-atom inter- 
action potential at small distances, see Eq. (j^. 

Ref. W4 We 

He(2 i5)-He(2 This work 3.912 7(5) 555.86(5) 

He(2 35)-He(2 35) Chen and Chung @ 3.3006 314.18^ 

This work 3.305 2(5) 314.44(5) 



■Multiplied by the factor | to correspond to the theory of Ref. 
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TABLE IX. For Re{2^S), values of -R^VAtD{R,e), where VA.tD(-R,e) is the atom-wall poten- 



tial, for values of e 


corresponding 


: to several types of dielectric, 


cols. 2-4, and in 


col. 5 values of 


-R^VAtMiR) for a 


perfectly conducting wall. The coefficient C3{2^S) is 2.67121. 




R 


e = 2.123 


e = 2.295 


e = 3.493 


e = oo 


10 


0.95339 


1.04221 


1.47123 


2.65990 


15 


0.95029 


1.03882 


1.46644 


2.65455 


20 


0.94739 


1.03564 


1.46194 


2.64938 


25 


0.94463 


1.03262 


1.45768 


2.64439 


30 


0.94200 


1.02975 


1.45361 


2.63956 


50 


0.93244 


1.01928 


1.43883 


2.62159 


70 


0.92395 


1.00999 


1.42570 


2.60532 


100 


0.91253 


0.99750 


1.40805 


2.58320 


150 


0.89577 


0.97916 


1.38214 


2.55042 


200 


0.88087 


0.96286 


1.35913 


2.52098 


250 


0.86726 


0.94797 


1.33812 


2.49381 


300 


0.85463 


0.93415 


1.31863 


2.46833 


500 


0.81082 


0.88624 


1.25109 


2.37768 


700 


0.77415 


0.84615 


1.19462 


2.29898 


1000 


0.72765 


0.79532 


1.12305 


2.19547 


1500 


0.66478 


0.72660 


1.02637 


2.04896 


2000 


0.61386 


0.67095 


0.94810 


1.92470 


2500 


0.57109 


0.62421 


0.88236 


1.81635 


3000 


0.53433 


0.58405 


0.82587 


1.72023 


5000 


0.42573 


0.46540 


0.65887 


1.41942 


7000 


0.35344 


0.38641 


0.54755 


1.20473 


10000 


0.28060 


0.30682 


0.43519 


0.97660 


15000 


0.20730 


0.22670 


0.32187 


0.73525 
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20000 0.16342 0.17873 0.25390 0.58540 

25000 0.13443 0.14702 0.20893 0.48434 

30000 0.11395 0.12463 0.17715 0.41205 

50000 0.07032 0.07692 0.10938 0.25591 

70000 0.05067 0.05543 0.07883 0.18478 

100000 0.03565 0.03899 0.05545 0.13013 
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TABLE X. For Re{2^S), values of -R^VAtD{R,e), where VAwiR,^) is the atom-wall poten- 



tial, for values of e 


corresponding 


: to several types of dielectric, 


cols. 2-4, and in 


col. 5 values of 


-R^VAtMiR) for a 


perfectly conducting wall. The coefficient C3{2^S) is 1.90092. 




R 


e = 2.123 


e = 2.295 


e = 3.493 


e = oo 


10 


0.67644 


0.73946 


1.04384 


1.88963 


15 


0.67336 


0.73608 


1.03907 


1.88428 


20 


0.67047 


0.73292 


1.03459 


1.87912 


25 


0.66773 


0.72992 


1.03036 


1.87413 


30 


0.66512 


0.72707 


1.02633 


1.86931 


50 


0.65566 


0.71671 


1.01169 


1.85142 


70 


0.64728 


0.70755 


0.99875 


1.83525 


100 


0.63606 


0.69527 


0.98140 


1.81333 


150 


0.61966 


0.67733 


0.95607 


1.78095 


200 


0.60516 


0.66146 


0.93368 


1.75197 


250 


0.59196 


0.64703 


0.91333 


1.72529 


300 


0.57977 


0.63370 


0.89454 


1.70030 


500 


0.53788 


0.58789 


0.83004 


1.61162 


700 


0.50336 


0.55016 


0.77695 


1.53484 


1000 


0.46046 


0.50328 


0.71103 


1.43436 


1500 


0.40433 


0.44195 


0.62481 


1.29413 


2000 


0.36074 


0.39433 


0.55783 


1.17825 


2500 


0.32560 


0.35593 


0.50380 


1.08034 


3000 


0.29654 


0.32419 


0.45909 


0.99640 


5000 


0.21739 


0.23769 


0.33709 


0.75412 


7000 


0.17046 


0.18640 


0.26458 


0.60138 


10000 


0.12784 


0.13981 


0.19860 


0.45715 


15000 


0.08946 


0.09785 


0.13907 


0.32318 
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20000 0.06848 0.07490 0.10649 0.24852 

25000 0.05536 0.06055 0.08610 0.20136 

30000 0.04641 0.05076 0.07219 0.16904 

50000 0.02810 0.03074 0.04372 0.10257 

70000 0.02012 0.02201 0.03131 0.07350 

100000 0.01411 0.01543 0.02195 0.05154 
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FIGURES 
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FIG. 1. Comparison of measured values (A, B) and the upper and lower bounds of Glover and 
Weinhold (shaded region) p^] with the present calculation of the static polar izability for He(2^5). 



The point A is the measurement of Crosby and Zorn [18| and the point B is that of Ekstrom et 



al. 1 16] determined in combination with measurements from Molof et al. |17|. 
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2 4 6 8 10 

3 

atom-wall distance (10 a.u.) 

FIG. 4. Potentials R^V{R) for Re{2^S) atom-wall interactions. The labels A, B, and C 
spond respectively, to dielectric constants e of 2.123, 2.295, and 3.493. 
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